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ABSTRACT 

Utilizing non-LTE time-dependent radiative-transfer calculations we investigate the impact 
of mixing and non-thermal processes associated with radioactive decay on SN Ilb/Ib/Ic light 
curves and spectra. Starting with short-period binary models of <J 5-Mq helium-rich stars, 
originally 18 and 25 M Q on the main-sequence, we produce 1.2 B ejecta which we artificially 
mix to alter the chemical stratification. While the total 56 Ni mass influences the light-curve 
peak, the spatial distribution of 56 Ni, controlled by mixing processes, impacts both the multi- 
band light curves and spectra. With enhanced mixing, our synthetic light curves start their 
post-breakout re-brightening phase earlier, follow a more gradual rise to peak, appear redder, 
and fade faster after peak due to enhanced 7-ray escape. Non-thermal electrons, crucial for 
the production of He I lines, deposit a large fraction of their energy as heat, and this fraction 
approaches 100% under fully ionized conditions. Because energy deposition is generally lo- 
cal well after the light-curve peak, the broad He I lines characteristic of maximum-light SN 
Ilb/Ib spectra require mixing that places 56 Ni and helium nuclei to within a 7-ray-mean-free- 
path. This requirement indicates that SNe lib and lb most likely arise from the explosion of 
stripped-envelope massive stars (main-sequence masses <J 25 M Q ) that have evolved through 
mass-transfer in a binary system, rather than from more massive single Wolf-Rayet stars. In 
contrast, the lack of He I lines in SNe Ic may result from a variety of causes: A genuine helium 
deficiency; strongly-asymmetric mixing; weak mixing; or a more massive, perhaps single, 
progenitor characterized by a larger oxygen-rich core. Helium deficiency is not a prerequisite 
for SNe Ic. Our models, subject to different mixing magnitudes, can produce a variety of SN 
types, including lib, lie, lb, and Ic. As it is poorly constrained by explosion models, mixing 
challenges our ability to infer the progenitor and explosion properties of SNe Ilb/Ib/Ic. 
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1 INTRODUCTION 

While the number of peculiar transients has grown in the last 
decade, much debate still surrounds the generic properties of stel- 
lar explosions, including core-collapse supernovae (CCSNe) and 
their progenitors. Ultimately, we aim to make an accurate predic- 
tion of the explosion energy, explosion morphology, ejecta mass, 
remnant mass, nucleosynthetic yields etc. of a massive star given 
its mass, rotation rate, metallicity on the main sequence, and bina- 
rity. At present, the exercise is done the opposite way - for any new 
SN a model is crafted to match observations with various levels of 
success. To reach a higher level of consistency, considerable work 
is still needed in stellar evolution, stellar explosion, and radiative- 
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transfer modeling of CCSNe. The present work is devoted to this 
later domain. 

There is strong evidence that SNe II-P are linked to red- 
supergiant star explosions, either through direct progenitor identi- 
fication on pre-explo s ion images or through light curve modelling 
dFalk & Arnetd|l977l: iLitvinova & Nadezhinl|l985l; lEastman et al.l 
1 1994 IUtrobinll2007l: iKasen & Wooslevl|2009t iDessart et al.ll2oTcT 
bessart & Hillied 12011 at iBersten et al.l t201lh . The failure to de- 
tect t he progenitor s of nearby SNe Ib/c on pre-explosion im- 
ages ISmartJ 120091) prevents such a robust identification. How- 
ever, light curve (|Ensman & WooslevFl 988; Shig evama et al.ll99Ct 
IDessart et al.ll2011 ) as well as spectral dSwartz et al.l 1993 ) mod- 
elling favors low-mass ejecta likely resulting from the explosion of 
essential ly hydrogen-less cores of moderate mass evo lved in binary 
systems jPodsiadlowski et al.ll992llYoon et al.l201fj) . The surpris- 
ingly high occurrence rates of SNe lb and Ic versus II-P in existing 



© 2010 RAS 



2 Luc Dessart et al. 



surveys dSmith et aljEoi l]), together with the non-detection of SNe 
Ib/c prog enitors so far, gives fu rther support to such binary-star pro- 
genitors (Eldridge et al.l2008l) . 

iGeorgv et al.l d 2009) have studied the evolution of single mas- 
sive stars and find a range of metallicities and initial masses for the 
production of helium-rich and helium-deficient Wolf-Rayet (WR) 
stars at death. Assuming a successful explosion follows the collapse 
of the core, they infer lb to Ic SN rates that corroborate the ob- 
servations as well as the expectations from binary-star evolutionQ 
However, this low-mass binary-star population is distinct from the 
recorded WR star population, which is dominated by higher-mass 
objectfl that shed their envelope through radiatively-driven winds 
rather than Roche-lobe overflow, and endowed with a higher mass 
at collapse. Such properties at death dCrowtherll2007l) do not seem 
comp atible with what makes the bulk of SNe Ib/c dDessart et al .1 
l201lh . 

Unstable 56 Ni nuclei, created by explosive nucleosynthesis, 
are ultimately responsible for the long-term brightness of CCSNe. 
In Type I CCSNe, the radioactive-decay energy from 56 Ni and its 
daughter nucleus 56 Co also powers the peak of the light curve, 
whose evolution is controlled by the dynamic radiative diffus ion 
of the heat released at depth by decay dColgate & McKedll969h . 

As the densities decline with age, the assumption that ra- 
dioactive decay energy is deposited only as heat becomes unten- 
able. Decaying nuclei produce 7-rays that Compton-scatter with 
free and bound electrons, producing a hig h-energy tail to the oth - 
erwise Maxwellian electron distribution dSpencer&Fanolll954l) . 
These non-thermal electrons are a source of heat, critical for the 
light curve, and a source of non-thermal excitation and ionization. 
The latter affect SN spectroscopic signatures, and are critical for the 
produ ction of the optical He I li nes that characterize the Type lb SN 
class dLucvll 1 99 ll ; fs wartj 1 99 ll) . They also play a considerable role 
in determin ing the excitation/ionization of the gas during the neb- 
ular phase dKozma & FranssorJll992l . Il998al lbl). Because of such 
compli cated physics, in addition to non-LTE and time-dependent 
issues dDessart &HidTieil 2008, 2010), estimating the exact compo- 
sition of these ejecta, and in particular the helium abundance, is a 
challenge. 

In the last decade, the observational divers i ty of SNe Ib/c 
has been documented dRichardson et alj 120061 ; iModiaj |2007| ; 
iDrout et al] |20 1 ll) . Several ne arby events received special atten- 



tion, such as the SN lb 200 8D dSoderberg et al.l2008l;lModiaz et al 



20091) or the SN Ic 19941 dFilippenko et al.lll995l;lRichmond etal 



199q ; IClocchiatti et al.l 1996 ). While helium is undoubtedly present 
in SNe lb, it is unclear whether the lack of obvi ous He I lines in 
SNe Ic spectra is evidence fo r helium deficiency dFilippenko et al .1 
ll995l : IClocchiatti et al.ll996T) . Line overlap and blanketing prevents 



1 In all these stellar-evolution studies, the definition of a SN lib, lb, and 
Ic is based on the progenitor surface composition which does not directly 
translate into a SN type due to the subtleties of atomic physics and ra- 
diative transfer. Much ambiguity thus surrounds the definition of a SN lb 
or Ic, which is distinct for SN observers, SN modelers, and modelers of 
stellar-evolution. However, with the exception of those evolved at twice- 
solar metallicity, single-star models have final masses >10Mq, and thus 
are much larger than those of pre-SN "hydrogen-less" cores produced from 
binary-star evolution. 

2 The reason is simply that WR stars are searched for using a 
luminosity/mass-loss criterion, reflected by the strength of their optical 
emission lines relative to the overlapping continuum. Only highly luminous 
WRs with dense stellar wind s, which are prim arily single massive stars, are 
detected in such surveys (see lCrowtheill2007l for a review). 



a definite resolution of this problem based on a single-epoch opti- 
cal or near-IR spectrum. The task instead requires a full modelling 
of the light curve and spectra, from photospheric to nebular phase, 
since the successful model must support observations at all times 
and wavelength ranges. Such results must also be confronted with 
the predictions from single- or binary-star evolution. 

So far, CCSN radiation modeling has been split into two dis- 
tinct efforts, one devoted to the light curve and the other, usually us- 
ing entirely different assumptions, to the spectra. Light-curve mod- 
elling of SNe Ib/c has been performed based on theoretical models 
of the hydrogen-less cores produced through bina ry-star evolution 
dEnsman & Wooslevll988l ; lshigevama et aljl990h . Because it is so 
critical to the problem, time-dependent radiation transport is treated 
but the gas is assu med in LTE and opacit ies are approximated. Us- 
ing this approach, llwamoto et al. I dl994h modelled the SN 19941 
light curve to infer a record-low ejecta mass of 0.9 Mq for a CCSN. 

In contrast, spectroscopic modelling typically allows for de- 
partures from LTE but assumes steady-state radiation transport 
and is limit e d to th e ph otospheric layers exc lusively. In this spirit, 
ISauer et al.l d2006l) and iTanaka et all d2009r) used a Monte Carlo 
technique to model multi-epoch spectra of SNe 19941 and 2008D, 
respectively. A similar photospheric steady-state approach assum- 
ing an homogeneous composition and a diffusive inner boundary 
(but with a l lowan ce for non-thermal effects) was performed by 
iBaron et al.l dl999l) on SN 19941. The assumption of steady-state 
makes the modelling more practical since one can adjust model 
parameters any number of ways in order to fit the spectrum but 
this limits both the accuracy and the consistency of the simulation. 
The neglect of time-dependent effects impacts the ionization and 
temperature structure; the prescribed diffusive inner-boundary flux, 
which becomes inadequate as the ejecta thins out (i.e., as early as 
the peak of the light curve), compromises the accuracy of synthetic 
spectra longward of the visual band. In addition, the piecemeal ap- 
proach of these studies, which compute separately light curves and 
spectra, is not globally consistent: the best-fit model to the observed 
light curve may conflict with multi-epoch spectra, and vice- versa. 

To remedy such shortcomings, we have developed a new non- 
LTE line-blanketed time-dependent radiative-transfer modelling 
approach that simultaneously computes the multi-band light curves 
and spectra with the capability to cover continuously from early 
times when the ejecta is optically thick to late tim es when it is thin 
dDessart & Hillie j|20 1 Ol ; iHillier & Dessartll2012l) . We have applied 
this method to study the early-time evolution of SN 1987A, and 
extended this work to SN e Il-Plateau up to 1000 d after explosion 
dDessart & Hillie J201 1 ah . 

More recently, to gauge the conditions under which one may 
see H I and He I lines in the optical and near -IR ranges, we stud - 
ied the early-time evolution of SNe Ilb/Ib/Ic dDessart et alj|201lh . 
These full-ejecta simulations were based on realistic binary-star 
evolution models dYoon et al.l201(il , Woosley et al. in prep), evolved 
until iron-core collapse, and exploded by means of a piston with al- 
lowance for explosive burning. These simulations indicate that as 
little as 0.001 Mq of hydrogen with surface mass fraction as low 
as 0.01 can produce a strong Ha line and a likely classification as 
lib (or lie as we argue in this paper). Contrary to myth, we also 
find He I lines can be produced in the absence of 56 Ni, but their 
presence is short lived and hinges on a very high helium mass frac- 
tion in the ejecta. From the general agreement of our synthetic light 
curves with observations, we argued that SNe Ilb/Ib/Ic are gener- 
ally low-mass ejecta, likely resulting from low-mass massive stars 
evolved in interacting binaries, as generally proposed. 

In this work, we repeat a sample of such SN Ilb/Ib/Ic simu- 
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lations but this time we include the treatment of non-therm al ef- 
fects a ssociated with 7-ray-energy deposition. As shown by iLucvl 
dl99lh . non-thermal effects are crucial for understanding the pro- 
duction of He I lines in SN lb spectra. However, with advances 
in computational power and the availability of better atomic data, 
we relax many of the assumptions made in his treatment, such as 
steady-state (applied to single-epoch spectra), not fully non-LTE, 
neglect of line-blanketing, and the use of fully mixed (i.e., homo- 
geneous) eje cta. Our paper is also a con tinuation of the earlier ex- 
plorations of IWooslev & Eastman! d 19971) on the role of mixing and 
non-thermal processes in SN lb and Ic ejecta. 

To investigate the effect of mixing on CCSN light curves and 
spectra originating from our binary-star progenitors, we use a "box- 
car" algorithm to modulate the species distribution in the ejecta and 
hence mimic m ulti-dimensional effects associated with massive- 
star explosions (kifonidis et alj|2003l . 120061 : Ijoggerst et alj|2009l ; 
lHammer et al.ll201oh . Our crude mixing prescription only affects 
the distribution of the elements — the density structure remains 
unchanged. While some mixing seems to occur in all SNe Ib/c, we 
find that it must be very efficient in helium-rich ejecta to yield a SN 
lb - moderate mixing, even in helium-rich ejecta, will still produce 
a SN Ic spectrum. This sets strong constraints on the nature of these 
events, both concerning the pre-SN star evolution and the explosion 
physics. 

This paper is structured as follows. In the next section, we 
discuss the various issues that surround mixing in CCSN ejecta 
(origin, nature, effects). In Section [3] we briefly summarize our 
modelling approach, focusing primarily on the treatment of non- 
thermal pro cesses since this aspect is an a d dition to the metho d 
described inlDessart & Hilliend2010l,l2011al) ; lDessart et al.ldioTlh ; 
Iffillier & Dessartj J2012I) " We then describe in Section [4] the ther- 
mal and non-thermal effects associated with 7-ray energy depo- 
sition and the production of Compton-scattered high-energy elec- 
trons, and how these affect the formation of He I lines (Section[5]l. 
Section [6] describes the ejecta properties, in particular the photo- 
spheric properties, and how these are reflected in the light curves. 
We then discuss in Section [7] the spectral properties of our grid 
of SN Ib/c models and their sensitivity to the efficiency of ejecta 
mixing. In Section[8] we compare these synthetic spectra to obser- 
vations in order to gauge the adequacy of our models. We finally 
discuss the relevance of these results and present our conclusions 
in Section|9] 



2 THE NEED FOR MIXING 

Before presenting detailed radiative-transfer models that include a 
treatment of non-thermal effects, we start by motivating the need 
for mixing, since it modulates the magnitude of these processes 
through the local variation of the 56 Ni mass fraction. 

The 7-ray mean free path £ 7 is defined by £ 7 = Xjn^p, where 
k 7 is the mass-absorption coefficient to 7-rays and p is the mass 
density. For a uniform ej ecta composition, k~, is a c onstant^ which 
we set to 0.03 cm 2 g _1 dKozma & Franssonll 19921) . so that l~, just 



3 At the energies of interest, gamma-rays lose most of their energy through 
Compton scattering. While the scattering is forward peaked, the gamma- 
ra ys loose little energy until they undergo a large angle scattering. As shown 
bv lSwar tz et al. 1 1995j) accurate answers for the gamma-ray transport can be 
obtained by treating gamma-ray energy deposition as a purely absorptive 
process. 



scales inversely with p. The density distribution in WR-star explo- 
sions is well descr ibed above ~3000km s -1 by a power law with 
exponent n ~ 8 fensman"& Woosley|[T988l ; iDessart et alj|201ll) 
and is thus a steep declining function of radius (or equivalently 
velocity). Even with the corresponding strong increase in i 7 with 
radius, the ejecta density is very large at early times so that l-y/R 
remains generally small. Hence, 7-rays cannot travel the distance 
that separates the inner ejecta layers, where the 56 Ni was explo- 
sively produced, and the optical photosphere where the bulk of the 
SN radiation escapes^ Such deeply seated 7-rays will completely 
"thermalize" and trigger the formation of the heat wave that will 
eventually affect the light curve. However, they will induce no vis- 
ible non-thermal effect at early times. 

In homologously-expanding ejecta, the density in any mass 
shell varies with the inverse cube of the time since explosion, and 
thus the optical depth to the surface of the ejecta varies as the in- 
verse square of the time. Hence, deeply emitted 7-rays will eventu- 
ally be able to directly influence the bulk of the ejecta mass. How- 
ever, this will happen no earlier tha n about 2 months af ter explo- 
sion in unmixed SN Ilb/Ib/Ic ejecta IDessart et aT1l201 ll) . and thus 
after the medium is thin to optical photons. Using as a reference 
the spherically-symmetric eje cta produced from t he explosion of 
the Bmil8mf4p41zl model of lDessart et al. I J201 lb, we illustrate in 
Fig.Q](left plot) the depth and time variation of l-y/R. This quantity 
is indeed less than unity below 10000 km s~ for up to ~30 d after 
explosion, so that the 7-ray influence is essentially local (there is 
no "transport") for up to a month. An alternative method of exam- 
ining gamma-ray transport is to examine the radial optical depth, 
integrated inward from the outermost grid point down to r 7 = 2/3 
(Fig.|T] right plot). This gives the location above which gamma-ray 
photons can "freely" escape. 

The critical observation that requires mixing and the early 
influence of non-thermal processes on the observed radiation in 
SNe Ilb/Ib is th e persistence of He I lines throughout the pho- 
tospheric phase (Modjaz 2007). To produce non-thermal excita- 
tion/ionization early-on, 56 Ni nuclei have to be mixed into regions 
where helium atoms are present, although some helium will be 
mixed inwards too. Quantitatively, the width of the He 1 5875A line 
at early times in the standard SN lb 2008D implies 56 Ni has to 
be present in ejecta shells traveling at ~ 10000 km s _1 (what this 
implies for the SN explosi on energy depends on the ejecta mass, 
which is debated; see, e.g. jTanaka et alj|2009h . When we discuss 
mixing to large velocities, we mean the large velocities commen- 
surate with the observed line profile (half) widths. What clearly 
emerges from our modeling is that helium richness alone is not 
sufficient to ensure a SN lb classification; the " right" mixin g has to 
occur to allow efficient non-thermal excitation (ILuc vlll99ll) . In this 
respect, SNe lb set a very stiff constraint on the explosion physics 
through this mixing efficiency, as well as on the composition and 
mass of the ejecta in which this mixing operates. 

Mixing of 56 Ni out to large radii/velocities in CCSN ejecta 
may take various forms but it is fundamentally associated with 
multi-dimensional effects. Mixing may arise from a large scale 
asymmetry of the explosion associated with a jet/bipol ar/unipolar 
morphology, as in the collapsar model dWooslevlll993l) , an intrin- 
sic large-scale asymmetry of the shock in an explosion through 



4 The lower the main sequence mass the smaller the (oxygen-rich) mass 
buffer between these two regions. Intuitively, we anticipate such 7-ray- 
transport to occur sooner in lower-mass ejecta from lower-mass main- 
sequence stars in binary systems. 
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Figure 1. Left: Evolution of the radial variation of the 7-ray mean free path, normalized to the local radius in the SN lib model Bmil8mf4p41zl of lDessart et alj 
1201 II) . A color coding is used to distinguish different epochs, selected at fixed intervals between 1 . 1 and 40 d after explosion. The dashed line refers to Z 7 = R, 
which is a representative threshold for the non-local influence of 7-rays. Right: Same as left, but now showing the 7-ray optical depth integrated inwards from 
the outer boundary of the ejecta. The dots indicate the location of the electron-scattering photosphere (i.e., r ea = 2/3). 



the neutrino mechanism dScheck et al . 2006; Marek & Jankal2009T) 
or the magneto-rotational mechanism dLeBlanc & Wilson! [l97(il: 



Bisnovatvi-Kog an et alj|197fj; iKhokhlov et alJI 199SI : iMaeda et al.1 
20021 : iBurrows et al.l l2007bllal iDessart et al.l 120081) . Mixing may 
also stem from the small-scale structure of the SN shock, as- 
sociated with convection and slowly -developing Rayleigh - Taylor 
and Kelvin-Helmho l tz instabilities iKifonidis et alj 12002 , l2006l : 
|joggerstetal]|2009l : iHammer et ail |2010| : iNordhaus et al.l l2010h . 
Note that even with mixing, the bulk of the 56 Ni mass remains 
deeply seated because mixing only affects the species distribution 
while the mass density is essentially unaffected and steeply rises 
towards the inner ejecta. 

Because the 7-ray mean free path tends to be much smaller 
than the scale of the system at early times (Fig [TJ, these different 
aspherical morphologies behave very differently for the purpose of 
non-local 7-ray transport and non-thermal excitation/ionization of 
the gas. To illustrate this, we use the model shown in Fig.Q]at an 
age of 26.8 d (which typically corresponds to the light curve peak 
in a SN lb) and study the geometrical properties imposed on mixing 
to generate non-thermal effects at the photosphere at that time. For 
this model, t 7 = 2/3, is located at 7700 km s _1 . 

The definition of the 7-ray mean free path Z 7 contains no in- 
formation on the radial or angular dependence of the escape prob- 
ability of 7-rays at a given time from a given location in the ejecta. 
A better quantity is the angle-dependent 7-ray optical depth, whose 
variation dr 7 over a step size ds « R p hot is equal to «: 7 pds (the 
photospheric radius i? p hot is used as a representative scale for the 
system). 

Using the tangent-ray meth od to compute optical-depth inte- 
grals along multiple directions jMihalaslll978l) . we compute the 
location of the r 7 = 2/3 surface as seen from a selection of 
heights/radii/velocities in the SN ejecta[f| This surface bounds the 
ejecta volume over which a 7-ray source may cause non-thermal 

5 This calculation uses the model density distribution, interpolated on a 
logarithmic radial grid, with 400 radial points to ensure a good resolution in 
particular for lateral directions. The tangent-ray method is fast since it does 
not require any interpolation; all ray quantities are calculated at grid points. 
Depending on the location of the 7-ray source, we adjust the location of the 
outermost grid point to be at 20000 or 40000 km s ~ 1 . 



effects. In Fig. [2] we show results for a 7-ray-emission site that is 
below, at, or above the 7-ray photospherelf] 

For locations below the 7-ray photosphere, the volume of in- 
fluence is essentially a sphere centered on the emission site. As the 
emission site is brought outward nearer the photosphere, this sphere 
gets stretched on its outer side, bulges somewhat laterally, while its 
inner side stays very close to the emission site. As the emission site 
crosses the photosphere, the r = 2/3 surface breaks open along the 
radial direction, reaching all the way to infinity (at this time, 7-rays 
would be observable along this beam). In comparison, the opening 
angle of that surface remains rather small, and wraps tightly to the 
emission site at depth. As the emission site is brought further out, 
the opening angle broadens, eventually reaching the opposite side 
of the ejecta for the outermost locations. Even in this situation, 7- 
rays do not penetrate significantly the ejecta layers interior to the 
emission site. 

We see that the influence of 7-rays is highly angle dependent, 
with a strong bias for escape in the radial direction for locations at 
and above the photosphere. This results from the contrast in scale 
height between the radial direction (i.e., ~ R/n for a power law 
density distribution with exponent n) and the lateral direction (i.e., 
~ 2R/^/n). At early times, Z 7 is so small that 7-rays will only 
influence their immediate vicinity. Thus to allow non-thermal ion- 
ization/excitation of the bulk of the He atoms (preferentially lo- 
cated in the outer ejecta), the 56 Ni and 56 Co 7-ray sources need 
to be thoroughly mixed with the material in the helium-rich shell. 
By thorough mixing, we mean that fingers or clumps of 56 Ni will 
need to be sufficiently densely distributed so that the inter-clump 
distance is on the order of l~ t . At a radial distance R, this translates 
into an inter-clump angular distance A9 ~ LJ R. For the example 
shown in Fig.[2]at 26.8 d after explosion, clumps with a thickness 
AR and spaced every ~10° would allow 7-rays to influence the 

6 This exercise is illustrative and thus the message conveyed is not altered 
by the somewhat arbitrary choice for r 7 . In practice, using 2/3 yields the 
volume over which ~50% of the 7-rays are absorbed, while non-thermal 
processes may be sufficiently strong with only 10% of such 7-rays. In that 
case, the edge of that volume would correspond to t 7 = 2.3 from the 
emission site and the contours shown in Fig. [2] would correspond to 14.5 d 
after explosion rather than 26.8 d. 
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Figure 2. Slices through a SN ejecta (containing its center of mass; we use model Bmil8mf4p41zl of iDessart et alj 1201 ll) at 26.8 d after explosion) and 
illustrating the area (shown in gray) that is within a 7-ray optical depth of 2/3, assuming various heights along the vertical axis for the 7-ray emission site 
(shown as a black dot). We assume the SN ejecta has a spherically-symmetric mass distribution, so that this area corresponds to a volume of revolution about 
the vertical axis. Its morphology changes drastically as the 7-ray-emission site crosses the 7-ray photosphere (dash-dotted line) located at 7700km s _1 at this 
time (where the optical photosphere is at that time depends on the model; see Section[6). For emission sites near the 7-ray photosphere, 7-rays are incapable 
of reaching regions of similar radius/velocity but at different latitudes/longitudes: the larger 7-ray mean free path in the radial direction makes it the preferred 
direction of escape. 



Table 1. Progenitor and ejecta properties for models Bmf4p41x[l^l] and Bmf5p09x[l^t]. In our nomenclature, "B" stands for binary (evolution), "mfXpYZ" 
for the final star mass X.YZ in Mq (i.e., at the time of explosion). The suffix xl-x4 refers to increasing levels of mixing. Since mixing affects the ejecta 
distribution but not the cumulative ejecta masses of species, the yields are the same within a mixing sequence. The remnant mass is 1,48 Mq (1.5 Mq) 
for model sequence Bmf4p41x[0^l] (Bmf5p09x[0^t]), and all model ejecta have the same kinetic energy of 1.2 B. For both series, the pre-SN evolution is 
computed at solar metallicity. In model sequence Bmf4p41x[0— 4], the inner 4 zones of the KEPLER ejecta are infalling. We trim these regions, and because 
of the different levels of mixing, this causes slightly different 56 Ni cumulative masses, in the range 0.157 to 0.177 Mq. Note that in the subsequent evolution 
with CMFGEN, the steep density variation in the innermost ejecta layers of model Bmf5p()9xn causes numerical diffusion issues, which are worse in cases of 
weaker mixing. 



Model 


Mi 


M t 




M H 


M He 


M C 


M N 


M 


M Ne 


M-Mg 


M Si 


Mp c 


M 56 


Ni 




[M Q ] 


[M Q ] 


[M Q ] 


[Mq] 


[M ] 


[Mq] 


[Mq] 


[Mq] 


[Mq] 


[Mq] 


[M Q ] 


[Mq] 


[M©] 


Bmf4p41x[0-4] 


18.0 


4.41 


2.91 


6.88(-3) 


1.63(0) 


1.63(-1) 


8.83(-3) 


5.36(-l) 


1.471-1) 


6.99(-2) 


1.05(-1) 


6.87(-3) 


1.57(-1) .. 


1.770-1) 


Bmf5p09x[(M] 


25.0 


5.09 


3.61 


1.39(-6) 


1.58(0) 


4.3K-1) 


8.56(-3) 


1.15(0) 


7.601-2) 


7.30(-2) 


7.01(-2) 


8.850-3) 


1.84(-1) .. 


2.20(-l) 



© 2010 RAS, MNRAS 000. MEl 



Luc Dessart et al. 



o.io 



Velocity [10 4 km s~'] 
0.31 0.41 0.59 0.78 



1.41 




Solid: ™Ni Bmf4p41xl 
Broken: He 

Dotted: 

/" 

/ 



2.0 2.5 3.0 3.5 4.0 

Lagrangian Mass [M ] 



0.03 



to 
S 



1.0 



Velocity [10 4 km s~'] 
0.22 0.40 0.54 0.70 



Solid: ™Ni Bmf5p09xl 
Broken: He 




2.0 3.0 4.0 

Lagrangian Mass [M ] 



1.32 



5.0 





1.0 



Solid: °°Ni Bmf4p41x4 
Broken: He 




1.0 



0.8 



0.4 



£ 0.2 



Solid: M Ni 


Bmf5p09x4 


Broken: He 




; Dotted: 

















2.0 



5.0 



2.5 3.0 3.5 4.0 2.0 3.0 4.0 

Lagrangian Mass [M s ] Lagrangian Mass [M s ] 

Figure 3. Helium (dashed line), Oxygen (dotted line), and 56 Ni (solid line; we show here the initial distribution unaffected by decay) mass distribution versus 
ejecta mass for Bmf4p41 (left) and Bmf5p09 (right) model sets. The suffix xl-x4 refers to increasing levels of mixing (see Tables[TJ{2}. In the top row panels, 
we also add the velocity information. 
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Table 2. Same a s for Table □ but now showing the surface composition for each model in each series. The models with suffix xO correspond to the unmixed 
ejecta studied in lDessart et alj feOllh and are repeated here for completeness and reference. We also quote the ejecta velocity that encompasses 99% of the 
56 Ni mass initially, i.e., V56n;, me P ea k bolometric luminosity I/ pea k ar >d me time to peak luminosity t p eak- 
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entire shell at R. Strongly asymmetric explosions propelling 56 Ni 
nuclei to large velocities in a "jet" or a strong bipolar/unipolar ex- 
plosion are not attractive alternatives since the non-thermal pro- 
cesses would only affect a restricted volume of the SN ejecta - 
primarily the high-velocity material associated with the "jet" that 
does not contribute much to the bulk of the SN radiation. 



3 NUMERICAL APPROACH AND SETUP 

The numerical procedure is the same as that use d in our recent 
simulations on SN 1987 A jDessart & Hilliej2010t). SNe II-Platea u 
jPessart & Hilliej|2011ah . and SNe Ilb/Ib/Ic jDessart et alj|201ll) . 
and will thus be brie fly summarized be l ow. A full description of the 
code is presented in iHillier & Dessarj d2012h . The first distinction 
with those works is that instead of adopting a unique ejecta com- 
position structure, we explore different stratifications produced by 
different levels of mixing. Secondly, we now treat explicitly the 
effects of non-thermal excitation and ionization — in our forme r 
studies the entire decay energy was treated as heat dLi et al.ll2012h . 



3.1 Pre-SN evolution and explosion models 

In this study , we e mploy a subset of the models discussed in 
Dessart et alj 1 2oTTh . whose pre-SN evolution was described by 
Yoon ct al. see also Woosley et al., in prep)0 These two 

models corresponds to their models 21 and 28. Namely, we focus 
on the binary-star models Bmil8mf4p41zl and Bmf25mi5p09zl, 
which result from the evolution of an 18 and an 25 M© star in 
a binary system at solar metallicity, with, at the time of explo- 
sion, a mass of 4.41 and 5.09 Mq, and a stellar radius of 12.3 
and 4.35 Rp, respectively. Th ese two models were exploded with 
KEPLER ( I Weaver et alj 1 1978b by driving a piston at the iron-core 
edge to yield an asymptotic ejecta kinetic energy of 1.2 B. Explo- 
sive burning was computed from first principles with an approxi- 



mate 13-isotope network that includes 56 Ni (for simplicity, 56 Ni 
is also the only unstab le isotope that we consider in this work). 
In iDessart et alj d201 ll) . the ejecta were unmixed, i.e., the com- 
position was taken directly from the results of stellar evolution as 



7 In future work, we will investigate a broader range of progenitor s, in- 
cluding single rotating and non-rotating WR s tars taeorev et al.l2012l) . and 
wider interacting binaries dClaevs et alj201 ll) . 



modified by the explosion. The radiative-transfer simulations were 
started with CMFGEN from such hydrodynamical inputs evolved 
to 1 d after explosion. To be more concise, we adjust our former 
model nomenclature and assign to these two ejecta the root-names 
Bmf4p41 and Bmf5p09. 

In this work, we wish to explore the effects of mixing and non- 
thermal processe s on SN radiation. W e selected these two models 
from the study of lDessart et al.l J201 ll) to investigate the role of non- 
thermal processes in helium-rich low-mass cores with and with- 
out hydrogen at their surface, thus covering potential SN types of 
Ilb/IIc and Ib/Ic. All numerical codes used here are ID, i.e., the stel- 
lar evolution is done in ID, as is the radiation-hydrodynamics of the 
explosion, and the subsequent detailed radiative-transfer modelling 
of the ejecta gas and radiation. So, to mimic the multi-dimensional 
effects that may be at the origin of mixing in CCSN ejecta (see 
Section[2j, we run a "boxcar" through the unmixed ejecta models 
Bmf4p41x0 and Bmf5p09x0, and increase the efficiency of mix- 
ing by increasing the mass range over which the material is mixed. 
This mixing is applied to the original model shortly after explosive 
burning ceases so that the evolution to 1 d properly accounts for the 
heating associated with this ejecta stratification (in particular for 
the 56 Ni heat source). We produce Bmf4p41 and Bmf5p09 models 
that cover unmixed (suffix xO), weakly mixed (xl and x2), mod- 
erately mixed (x3) and strongly mixed (x4) configurations. Mixing 
does not change the cumulative yields of models Bmf4p41x0 and 
Bmf5p09x0, shown in Table [T] but it gradually modifies the chem- 
ical stratification, shown for each model in Fig. [3] Table [2] gives 
additional information on the surface composition for each model(f| 

Mixing softens composition gradients. In the most mixed 
models, the ejecta is nearly homogeneous — any species abun- 
dant at one depth is present at all depths. For stars of increasing 
mass, a greater mass buffer exists between the ejecta base and the 
helium-rich outer shells. Progressing outwards from the innermost 
ejecta layers, one encounters an oxygen-rich helium-deficient shell 
(here bounded between 1.8 and 2.5 Mq), followed by a shell with 
a composition shared between He, C and O. In our models, only 
the outermost layers (<0.5Mq below the surface) are nearly pure 



8 Since the mixing is performed over a specified mass range, it can affect 
all the way to the surface so a more desirable mixing approach would be 
to use a velocity criterion when setting the range of mass shells that are 
homogenized. 
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helium (with a mass fraction that depends on the level of mixing). 
In model Bmf4p41x0-3, hydrogen is also present with a mass frac- 
tion in excess of 0.01, which is su fficient to produce Ha emission 
at early times dDessart et al.ll20lH) . 

The velocity of the ejecta shell that bounds 99% of the to- 
tal 56 Ni mass is a good indicator of the magnitude of mixing 
in our models. This velocity varies from ~ 1500 km s _1 up to 
<8000km s _1 (Table EJl. Note that the bulk of the mass of such 
unstable isotopes is still located at small velocities because the mix- 
ing, which affects the profile of species' mass fractions, leaves the 
steeply-rising density distribution at depth unaffected. 

3.2 non-LTE time-dependent radiative-transfer including 
non-thermal processes 

The critical improvement in this work over all past studies of SNe 
Ibc explosions is the simultaneous computation of the SN multi- 
color light curve and spectra with allowance for line-blanketing, 
non-LTE, time-dependence in the statistical, energy, and radiative- 
transfer equations, and non-thermal processes associated with high- 
energy electrons Compton scattered by the 7-rays arising from ra- 
dioactive d ecay. All ingredie n ts but the las t have been de s cribed 
in d etail inlDessart & Hilliej fcOloL l2011ab ; iDessart et all d201lh 
and iHillier & Dessartj d2012h . In these former works and in the 
present study, time dependence is considered both for the treat- 
ment of the radiation field (i.e., DJ/Dt and DH/Dt terms in the 
0th and 1st moment of the radiative-transfer equation) and of the 
gas (De/Dt and Dn/Dt terms appearing in the energ y and the 
statistical equilibrium equations; IDessart & Hillierll2008l ; we thus 
naturally account for time-dependent ionization). The new addi- 
tion is the treatment of non-t hermal processes throu gh a solution 
of the Spencer-Fan o equation I Spencer & Fanol 19541) using the ap- 
proach presented in lKozma & Franssonld 19921) . Assuming a 7-ray- 
energy deposition profile, we inject the corresponding energy in the 
form of electrons at 500 eV and compute the local degradation of 
such electrons as they collide with electrons and ions/atoms in the 
plasma. Once this degradation function is known, we can determine 
how this deposited energy is channelled into the various forms of 
heat, ionization and excitation. The non-thermal terms associated 
with ionization and excitation are entered as rates in the statistical- 
equilibrium equations and are thus fully integrated in the CMFGEN 
solution to the non-LTE time-dependent ionization and level popu- 
lations. A full description of the technique and an in-depth analysis 
of non-therma l effects in cor e-collapse SNe is presented in a com- 
panion paper dLi et al 1 l2012h . Here, we will limit the discussion of 
non-thermal effects relevant to understanding the formation of He I 
lines (Section[4]l in SNe Ibc. 

In practice, the non-thermal solver in CMFGEN is only used if 
the mass fraction of unstable isotopes at the photosphere is indeed 
sizeable, typically above 0.0001. This was decided for each model 
by switching on the non-thermal processes at selected times and 
checking if this produced any effect on the gas properties and/or on 
the SN radiation. As discussed in Section|2] / 7 in the 56 Ni-rich re- 
gions is small at early times in all but the most mixed models (suffix 
x4), so we switch on non-thermal processes at the start of the time 
sequence only in the x4 series. To save space and since the u nmixed 
models have already been discussed in lDessart et al. In, we do 
not show results for models Bmf4p41x0 and Bmf5p09x0; models 
Bmf4p41xl and Bmf5p09xl behave very similarly to them and are 
thus used as substitutes. 

For the 7-ray-energy deposition, we can either assume it is 
local or we can compute its distribution in the ejecta using a 7- 



ray M onte-Carlo transport code we developed dHillier & Dessarj 
l2012h . We find that the 7-ray-energy deposition is essentially local 
up to about 30 d after the light curve peak in all models. Because we 
focus on the photospheric-phase evolution all simulations shown 
here assume a local energy deposition. Non-local energy deposition 
is not crucial for the non-thermal effects discussed here since it 
merely stretches the deposition profile outward, with an enhanced 
escape at large velocity. It does, however, change the amount of 
energy trapped in the ejecta and thus the luminosity decline rate 
during the nebular phase. 

The model atoms a dopted for all simulations in this work are 
identical to that used in IDessart & Hillierl d2011ah . The complete- 
ness of the model atoms influences the spectral appearance but it 
does not alter the He I-line strength or the non-thermal effects we 
observe. The main impact of increasing the size of the model atoms 
in our simulations, in particular for metals, is to enhance the mag- 
nitude of line blanketing, which tends to make the spectral-energy 
distribution (SED) redder. This effect is generally weaker than ob- 
tained here through variations in composition caused by mixing 
(i.e., Ni, C/O versus He etc.), although it noticeably alters the col- 
ors after the light-curve peak. Since atomic data represents an im- 
portant part component of the work we do, we wish to provide the 
sources and references for the data we use. The sources of atomic 
data are varied, and in many cases multiple data sets for a given ion 
are available. In some cases these multiple data sets represent an 
evolution in data quality and/or quantity, while in other cases they 
represent different sources and/or computational methods. Com- 
parisons of models calculated with different data sets and atomic 
models potentially provide insights into the sensitivity of model re- 
sults to the adopted atomic atoms and models (although such calcu- 
lations have yet to be undertaken for SN). Oscillator strengths for 
CNO eleme nts were originally taken from iNussbaumer & Storevl 
dl98j| 19841) . These authors also provide transition probabilities to 
state s in the ion continuum . The largest source of oscillator data is 
from lKuruczl d2009. 2010); its principal advantage over many other 
sources (e.g., Opacity Project) is that LS coupling is not assumed. 
More recently, non-LS o scillator strengths hav e become available 
through the Iron Project dHummer et al.lll993l), and work d one by 
the atomic-data group at Ohio State University dNahar|2010l). Other 
impor t ant sou rce s of radiative data for Fe include Becker & Bulled 
dl992l.ll995al lbl). lNahaJ dl995l) . Energy levels have generally been 
obtained from NIST. Collisional data is sparse, particularly for 
states far from the ground state. The principal source for collisional 
da ta among low lyi ng states for a variety o f species is the tabulatio n 
bvlMendozal dl983h; ot her sources includ e lBerrington et al.l(fl985l). 
iLennon et all 1 1985b jLennon & Burke] dl994l) IShine & Linskvl 
dl974l) . lTava1dl997al lbl), Zhang & Pradhan dl995d ; ll995bl ; ll997b\ 
Photo ionization data is ta ken from the Opaci ty Project dSeatonl 
Il987b and the Iron Project dHummer et al.ll993l) . Unfortunately Ni 
and Co photoionization data is generally unavailable, and we have 
utilized crude approxi mations. Charge exchange c ross-sections are 
from the tabulation bv lKingdon & Ferl and ( 1996). 

Finally, we have identified a numerical-diffusion problem 
which arises from the repeated remappings performed at each step 
in a time sequence. These remappings are necessary to adjust the 
grid onto the optical depth scale, which changes due to expansion 
and recombination/re-ionization. Unfortunately, the ID hydrody- 
namical inputs we employ have strong composition and density 
gradients which are smoothed by such remappings. Since the model 
at a new time step uses the model at the previous time step, these 
remappings produce systematic shifts which become continually 
worse as we evolve the models to later time. This alters the total 



© 2010 RAS, MNRAS OOO.mfSTl 



Radiative transfer ofSN Ib/Ic ejecta 9 



mass of some species by 10%, but by as much as 30% for 56 Ni in 
unmixed models because Ni is most abundant in a narrow shell 
at the ejecta base where the density steeply rises. While this does 
not impact the discussion on non-thermal effects, it induces a slight 
shift between model light curves (most visible at light curve peak). 
In future modeling this problem, at least for models where the ve- 
locity of a given mass shell is constant, will be rectified by interpo- 
lating the density structures from the original hydrodynamic model 
grid. Radiation-hydrodynamics codes do not suffer from this short- 
coming because they use a high-resolution Lagrangian mass grid, 
but it is then the poor resolution of the optical-depth scale that af- 
fects their results. 



4 THERMAL AND NON-THERMAL EFFECTS 

ASSOCIATED WITH 7 -ray ENERGY DEPOSITION 

Before presenting light curves and synthetic spectra obtained for all 
model sequences, we first describe the effects associated with 7-ray 
energy deposition in two representative models - one with strong 
mixing (Bmf4p41x4) and one with weak mixing (Bmf4p41xl). 
These two models are characterised by the same total mass, the 
same cumulative yields (i.e., same total mass of individual species), 
the same density structure, and the same explosion energy (or ex- 
pansion rate). Hence both models have, at a given post-explosion 
time, the same mass-density structure (versus radius or velocity). 
The differences between them arise directly from differences in 
chemical stratification and in particular how the different 56 Ni 
distributions influence the heating/ionization/excitation of the gas. 
Further, mixing impacts the opacity and emissivity of the gas 
through changes in the chemical abundances. 

Using the highly mixed model Bmf4p41x4, we show in Fig. [4] 
the ionization fraction (x e ; black), the fractional energies going 
into heat SEh (blue), ionization SEi (green), and excitation SE e 
(red), and the normalised energy deposition rate (e 7 . n ; orange) ver- 
sus ejecta velocity at 10.3 (top), 29.3 (middle), and 63.0 d (bottom) 
after explosion. To gauge the potential effects of non-thermal elec- 
trons on the continuum and line formation regions, which are lo- 
cated at and above the optical photosphere, we show this location 
in each panel with a black dot. We find that at all three epochs, most 
of the energy originally in non-thermal electrons gets degraded into 
heat. Hence, even if injected non-fhermally, most of this energy 
ends up thermal here and directly influences the light curve. Near 
the photosphere, 10% of the total deposited energy goes into ion- 
ization and about half that goes into excitation. Because helium is 
the most abundant species throughout most of the ejecta, it is He I 
that benefits the most from this non-thermal ionization and excita- 
tion (dashed curves). 

In a complicated fully-coupled non-LTE situation, there is no 
strict separation between these various heat, ionization, and exci- 
tation contributions. For example, heat can indirectly cause photo- 
ionization from a thermal blackbody and collisional ionization with 
thermal electrons at large optical depth where the material is close 
to LTE. Furthermore, over-excitation of upper levels in an atom/ion 
will facilitate photo-ionization and thus indirectly influence the ion- 
ization state of the gas. These complications justify a full non-LTE 
treatment so that all important rates are allowed for in the determi- 
nation of the level populations, the electron density, and the gas 
temperature. This treatment must also allow for time-dependent 
terms in the statistical equilibrium equatio ns since these are know n 
to influence the ionization state of the gas l lDessart & Hillierl2008l) . 
and thus the fate of non-thermal electron energy. 
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Figure 4. Illustration of the fractional energies going into heat SE^ (blue), 
ionization SEi (green), and excitation 5E e (red) as a function of ejecta 
velocity and shown at 10.3 (top), 29.3 (middle), and 63.0 d (bottom) after 
explosion in model Bmf4p41x4. The fractional energies shown with solid 
lines denote the total contributions in each energy channel, while the cor- 
responding dashed lines give the ionization and excitation associated with 
He I only. We also show the ejecta ionization fraction x e (black), which we 
take here as the ratio of the total ion population to the combined total atom 
and ion populations, and the normalised decay-energy deposition profile at 



the time (e~ n , in units of erg cm 



" 1 , and normalized to its maximum 



value; orange). In each plot, the black dot represents the location where the 
inward-integrated Rosseland-mean optical depth is 2/3. 
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Figure 5. Same as Fig. \4\ but now for the very weakly mixed model 
Bmf'4p41xl at 29.4 d after explosion. The relative fractions of the energy 
deposited by radioactive decay are comparable to those obtained for model 
Bmf4p41x4 at the same time (middle panel of Fig. [4), although the lower 
ejecta ionization enhances somewhat the contributions to ionization and ex- 
citation (which primarily go to He I). However, the magnitude of these non- 
thermal processes is completely dwarfed at and above the photosphere by 
the very steeply declining energy deposition profile (orange). In fact, most 
of the decay energy is deposited in regions that are ionized, located below 
the photosphere (black dot), so that the bulk of the decay energy goes into 
heat in this unmixed model. 

The magnitude of the effects associated with each channel 
scales with the actual energy deposited by radioactive decay. In the 
highly mixed model (Fig.O, the decay energy is deposited not only 
in the inner ejecta but also in high-velocity regions at and beyond 
the photosphere. This extra thermal and ionization energy causes 
the SN to brighten and peak earlier compared to unmixed models 
(see Section|6](. As we shall see in Section|7] it also leads to a siz- 
able strengthening of He I lines in this model series. 

In the weakly mixed model Bmf4p41xl (shown at 29.4 d 
only), the actual channelling of non-thermal energy into heat, ion- 
ization, and excitation, is very similar (Fig.|5)- The magnitude and 
distribution of the energy deposited is, however, drastically differ- 
ent. It is negligible at, and above, the photosphere, since the bulk 
of the 56 Ni is confined to the inner ejecta layers located below 
2800 km s _1 . Since models Bmf4p41xl and Bmf4p41x4 have the 
same amount of 56 Ni at a given post-explosion time, the tempera- 
ture in those inner layers is much higher than in the contempora- 
neous model Bmf4p41x4. As a consequence, the bulk of the decay 
energy is deposited below the photosphere, in deep regions that at 
this time are ionized. Consequently, the bulk of the injected non- 
thermal energy not only ends up as thermal energy, but it also fails 
to influence non-thermally the line and continuum formation re- 
gions. 

Non-thermal ionization and excitation are thus favored when 
the decay energy is deposited in regions of low ionization, which 
also tend to be external to the photosphere. Mixing is thus critical to 
propel 56 Ni at large velocities (and in all directions; see Section|2j, 
in regions that on a short time scale will be above the photosphere. 
In unmixed ejecta, heating from the 56 Ni at depth causes a strong 
and delayed heat wave that causes the photosphere to be external 
to the deep-seated 56 Ni-rich regions jDessart et al.ll20lH) . As an 
aside, the dependency of non-thermal effects to ionization suggests 
that they may not be paramount in SNe la ejecta since the large 
56 Ni mass causes a large ejecta ionization at all times. 



Interestingly, we find that the strong mixing adopted in models 
Bmf4p41x4 and Bmf5p09x4 actually quenches the impact of non- 
thermal ionization/excitation on He I in the outer ejecta. Indeed, 
energy deposited at large velocities causes heating of these layers, 
which over time and th rough the ionization fre eze-out caused by 
time-dependent effects dPessart & Hill ier 2008), cause helium to 
be fully ionized in those outer regions (the ion that gets the largest 
share of non-thermal ionization and excitation energy at large ve- 
locities is indeed Hell in Fig.[4]l. In this context, a lower level of 
mixing with allowance for non-local energy deposition could have 
allowed helium to recombine, and thus a stronger non-thermal ex- 
citation of He I atoms. In the future, we will need to investigate 
ways of mimicking as close as possible the mixing properties of 
multi-dimensional explosion simulations. 



5 PROCESSES CONTROLLING THE He I LEVEL 
POPULATIONS 

In SN ejecta, departures from LT E are usually very strong at, and 
above, the photosphere (see, e.g. JPessart & Hillierl J201 lah for an 
illustration of the departure coefficients of H I, O I, and Fe II in a SN 
II model). Here, non-thermal effects drive the material even further 
from LTE. 

In this section, we describe in more detail how the He I level 
populations are controlled by showing the dominant rates in and out 
of the lower nine levels (Fig. [6}. To assist with the interpretation, 
we give a succinct summary of the properties of these levels and 
the associated bound-bound transitions in Table [3] 

Rates controlling the level populations fall into different cat- 
egories but are fundamentally collisional or radiative. The former 
include thermal and non-thermal (collisional) ionization and excita- 
tion, and the latter photoionization and recombination (bound-free/ 
free-bound transition), and radiative excitation and de-excitation 
(bound-bound). In each panel of Fig.|6]we display a given rate if it 
exceeds 5% of the total rate at any given depth in the model (see 
label and color coding). The rate associated with a line transition 
at frequency v u i corresponds to the quantity n u A u i(l — J Vul /S u i) 
where n u is the population of the upper level, A u i the Einstein 
coefficient, J Vul the mean intensity at the frequency v u i, and S u i 
is the line source function. This rate may be positive (implying a 
net decay of electrons from the upper level) or negative (implying 
a pumping of electrons from the lower state into the upper state) 
depending on how J Vul compares to S u i- 

Below ~ 4000 km s _1 recombination and photoionization 
rates appear to dominate. However at these depth we are almost 
in LTE and many processes contribute to setting the level popula- 
tions. The line rates are also large, but since we plot the net rates 
they have been analytically cancelled. Above ~4000km s _1 indi- 
vidual processes are (generally) no longer in detailed balance, and 
the level populations are controlled by a variety of processes. 

In our treatment of He I non-thermal ionization is only al- 
lowed from the ground state of He I and thus only directly affects 
the population of level Is 2 1 S (top-left panel) and the ground state 
of He II. This is a valid approximation since the ground state is the 
most populated. Unlike recombination and photoionization rates 
that tend to cancel each other, there is no significant reverse pro- 
cess, and so, despite the modest rate, the impact on He ionization 
is important. 

While many different non-thermal excitation routes are con- 
sidered, the most important is from the ground state (Is 2 X S) to 
the Is 2p 1 P° level. When the bound-bound transition to the ground 
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Figure 6. Illustration of the dominant fractional rates in (positive) and out (negative) for the first 9 levels of He I, ordered from left to right and top to bottom 
Is ls 1 S, ls2s 3 S, ls2s 1 S, ls2p 3 P°, lslp 1 ? , ls3s 1 S, lsSp 1 ? , and ls3d 3 D (we omit level ls3s 3 S). For this illustration, we use model Bmf4p41x4 
at light curve peak. For better visibility, we display the individual rates as a fraction of the total rates for a given level at a given depth. [See Section [6] for 
discussion.] 
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Figure 7. Bolometric light curve for the Bmf4p41 series shown with respect 
to the time of peak bolometric luminosity and assuming full 7-ray trapping. 
In pr actice, 7-ray escape starts a few weeks after the peak in such progen- 
itors l lDessart et alj|201ll) and is slightly greater for higher mixing magni- 
tudes (at 100 d, we find that the fraction of 7-rays that escape increases 
from ~25% to ~40% from Bmf4p41xl to Bmf4p41x4). [See Section[6]for 
discussion.] 

state is sufficiently optically thick, the ls2p 1 P° level decays to 
Is 2s X S producing He 1 20581 A. Other 1 P° levels are also signifi- 
cantly influenced by non-therm al excitations from the ground state 

- in the Bethe approximation (I van Regemorteil [ l962f) the rate of 
excitation to these levels is proportional to the oscillator strength. 

In contrast, level ls2p 3 P° is primarily populated through 
downward transitions from higher levels and thus the 10830 A 
line from transition He I(ls2p 3 P°-ls 2s 3 S) stems primarily from 
recombination following photoionization and non-thermal ioniza- 
tion. He I 10830 A and 20581 A dominate the cascade into levels 
2 and 3, i.e., ls2s 3 S and ls2s 1 S (red curves in the 2nd and 3rd 
panels from top-left). 

The distinction between the ls2p 1 P° and the ls2p 3 P° pri- 
marily arises because of the large difference in collision cross- 
section. The collision strength between levels with permitted tran- 
sitions scales as logi? (energy of electron) at high energies (i.e., 
E ^> ionization energy), while the collision strength of other tran- 
sitions is almost constant (transitions to Is 2s 1 S) or declines with 
E (transitions to Is 2s 3 P°)0 Two other factors also come into play 

— for He I the rate of non-thermal ionizations is somewhat larger 
than the rate of non-thermal excitations, and approximately 3/4 of 
the recombinations (excluding the ground state) are to triplet levels. 

An important issue to remember when discussing the influ- 
ence of various processes is that the size of a particular process 
does not necessarily give an indication of its importance. An obvi- 
ous example is bound-free processes for the ground state of He I. 
While the component rates for the process may be larger than other 
rates, they may be negligible in controlling the helium ionization 
structure when the He I ground-state continuum is in detail balance, 
and the two rates (photoionization and recombination) cancel with 
a high degree of accuracy. 

In the present case non-thermal excitation and ionizations lead 
(directly or via recombination and cascades) to an over population 

9 We utilized the online Los Alamos Atomic Physics Codes at 
http://aphysics2.lanl.gov/tempweb/lanl/ which are based on Hartree-Fock 
method of R. D. Cowan, to compute collision strengths for both H I and 
He I. 



Table 3. Selection of properties of He I lower levels and important transi- 
tions in between them. We present an illustration of the rates controlling the 
population of the lower levels of He t in Fig.|6]in model Bmf4p41x4 at light 
curve peak. We first give the level Id. number, the level name, its statistical 
weight g, its excitation energy E{ (in cm -1 ) above the ground state. We 
then list the main transitions between these levels that fall in the optical and 
near-IR ranges, giving the oscillator strength and the wavelength for each. 



Id. 


Level 


3 


Ei (cm -1 ) 


1 


Is Is !S 


1 


0.0 


2 


ls2s 3 S 


3 


159856.1 


3 


ls2s 1 S 


1 


166277.5 


4 


ls2p 3 P° 


9 


169087.0 


5 


ls2p 1 P° 


3 


171135.0 


6 


ls3s 3 S 


3 


183236.9 


7 


lsSs^ 


1 


184864.9 


8 


ls3p 3 P° 


9 


185564.7 


9 


ls3d 3 D 


15 


186101.7 


10 


lsSd^ 


5 


186105.07 


11 


IsSp 1 ? 


3 


186209.47 


Transition 


gf 


Wavelength [A] 




4-2 


5.399E-01 


10830.17 




8-2 


6.406E-02 


3888.64 




5-3 


3.294E-01 


20581.28 




11-3 


1.558E-01 


5015.68 




6-4 


6.953E-02 


7065.25 




9-4 


6.097E-01 


5875.66 




7-5 


4.703E-02 


7281.35 




10-5 


7.094E-01 


6678.15 




11-3 


1.558E-01 


5015.68 





of the n=2 levels (primarily the 2s states) in He I. The population of 
these state is set by the leakage back to the ground state which com- 
pensate for the non-thermal excitations and ionizations. When this 
occurs photo-ionizations from these states may become the domi- 
nant ionization process, even though this process would be negli- 
gible i n the absence of non-thermal processes. As noted by iLucvl 
dl99lh this can amplify the influence of non-thermal processes. 

Non-thermal processes act directly only on a few levels in 
He I. However non-thermal ionization affects all levels indirectly 
through its influence on the He ionization and hence on recombina- 
tion. Non-thermal processes are important throughout the evolution 
of SNe Ibc, including the nebular phase. 



6 LIGHT CURVE AND PHOTOSPHERIC PROPERTIES 

In our non-LTE time-dependent radiative-transfer simulations, we 
compute the emergent spectrum from a few A to a few hundred 
microns. The frequency integral of this emergent spectrum at each 
epoch yields the bolometric luminosity of each model. We show 
the evolution for the bolometric luminosity for the Bmf4p41x[l-4] 
model series in Fig. [7] using the peak of the light curve as the time 
origin (the Bmf5p09x[l-4] series behaves similarly and is thus not 
shown). Full 7-ray trapping is assumed at all times. 

Enhanced mixing causes an earlier and more gradual rise to 
the peak. The post-breakout plateau that characterises the unm ixed 
or weakly mixed models, as described in lDessart et al.ldioilh van- 
ishes for the hi ghest levels of mixing (this i s similar to what was 
found by, e.g.. Istiigevama & Nomotolll990l for SN 1987A). En- 
hanced mixing causes an earlier peak (merely a few days in this 
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series but as much as ~ 15 d in the Bmf5p09x4 model), and an ear- 
lier transition to 7-ray escape, making th e nebular luminosity drop 
earlier and faster (not shown here, but see lDessart et al .1201 lh . This 
effect is however weak because, irrespective of mixing, the bulk of 
7-ray energy deposition reflects the 56 Ni distribution that is cou- 
pled with the density, and which steeply increases towards the in- 
ner ejecta. As we discussed earlier in the paper, some numerical 
diffusion affects somewhat the total mass of unstable isotopes, in 
particular in unmixed models. Since the main focus of this paper is 
on the non-thermal effects on SNe Ib/c spectra, we defer to a forth- 
coming study a more quantitative description of mixing on SNe 
Ib/c light curves (Dessart et al, in preparation). 

The light curve behavior reflects the evolution of the photo- 
spheric properties (Fig. [8} and nearly exclusively the radius and the 
temperature. For both model series Bmf4p41 and Bmf5p09, mixing 
causes the same effect on the photosphere behavior and primarily 
reflects the induced change on the 56 Ni distribution. With enhanced 
mixing, the heating of the outer ejecta layers is increased and the 
diffusion time of that energy to the outer layers is shorter. Con- 
sequently, the photosphere migrates outward in radius earlier with 
enhanced mixing. The extra heat produces a higher temperature ini- 
tially, which thus causes the luminosity to rise earlier, even avoid- 
ing a post-breakout plateau for the highest level of mixing here. The 
reduced mass fraction of 56 Ni at depth causes a smaller tempera- 
ture rise associated with the heat wave. With increased mixing, the 
non-monotonicity of the photosphere trajectory in mass/velocity is 
reduced and even disappears for the highest level of mixing. In this 
case, the temperature, the velocity and the mass above the photo- 
sphere continuously and gradually decrease with time. The extra 
heat generated from the mixed radioactive material keeps the pho- 
tosphere further out in the ejecta, at larger velocities. Without mix- 
ing, the photosphere recedes in mass very fast and abruptly changes 
its properties when it suddenly feels the heat wave from decay. 

These effects are all fundamentally of non-thermal origin, but 
as described in the previous section, the energy from non-thermal 
electrons is always channeled primarily in the form of heat. In the 
low-ionized conditions above the photosphere, it represents ~70% 
of the local energy deposition, reaches 95-100% far above the pho- 
tosphere where the conditions are ionized (which results primarily 
from ionization freeze-o ut associated with time-dependent effects; 
iDessart & Hilliej|2008h . Below the photosphere, the heat chan- 
nel represents ~95% of the total energy deposited locally, but de- 
creases with time at and beyond the peak of the LC, i.e., ~92% at 
29.3 d and ~82% at 63.0 d after explosion. Since the energy deposi- 
tion occurs primarily at depth, the bulk of the non-thermal energy is 
deposited as heat. Decay energy thus influences primarily the light 
curve, allowing such Type I CCSNe to be bright a few weeks after 
explosion. 



7 INFLUENCE OF MIXING ON SYNTHETIC SPECTRA 

Variations in mixing cause changes in composition (Fig. [3}, light 
curve evolution (Fig. [7}, and in non-thermal processes (Figs.[4}j5]l. 
In Figs.[9]and[T0] we show the spectral evolution over the first 100 d 
for models Bmf4p41xl-4 and Bmf5p09xl-4, respectively (black 
curves). We overplot in each panel the corresponding spectra when 
He I lines are excluded from the formal solution of the transfer 
equation. All simulations are done with large, but not huge, model 
atoms (Section[3j note that these simulations represent nonetheless 
the most CPU-intensive non-LTE line blanketed simulations ever 
performed for SNe Ib/c), so the spectral colors and morphology are 



not fully converged. For the present qualitative presentation, this 
is sufficient but we will explore the dependency of our results on 
the adopted model atom in a forthcoming paper. We obtain three 
different types of spectral evolution, corresponding to different SN 
classification: 

• The weakly and moderately mixed models (i.e., Bmf4p41xl- 
2 and Bmf5p09xl-3) have a similar spec tral appearance as the cor- 
responding unmixed models presented in lDessart et al. I d201ll) . De- 
spite variations in mixing, which impact both 56 Ni and all other 
species, the differences in this set are subtle. Most important, none 
of these models show optical He I lines after about 10 days past ex- 
plosion, and when such He I lines are present early-o n, they arise 
exclusively from non-LTE effects jDessart et al.l201lh . In the opti- 
cal, model Bmf4p41xl (and Bmf4p41x2, not shown) shows strong 
Ha initially, which quickly vanishes to leave a spectrum without 
H I, He I lines, nor the strong Si II absorption at 6300A. In a strict 
sense, this spectral evolution would correspond to a SN lie, a SN 
subtype unobserved today. All other models in this set would be 
of Type Ic. We note that the presence of He I and Ha in these 
models at early times is influenced by the ionization freeze-out in 
the outer ejecta, a time-dependent effect jDessart & Hilliej|2008l) 
that dominates in th e absence of non-thermal excitation/ionization 
jDessart et alfeoill) . 

• In the moderately mixed model Bmf4p41x3 (middle panel), 
Hel5875A is visible but it is extremely weak and very narrow, 
disqualified for a standard SN lb spectrum. Nonetheless, we can 
clearly see that the increasing mixing starts to have a visible ef- 
fect on the spectra. Non-local energy deposition at late times com- 
pensates for the weak mixing and causes the non-thermal excita- 
tion/ionization of He 1 10830 A; this effect is however too weak to 
produce unambiguous He I optical lines. Because of the presence 
of Ha at early times, this spectral evolution corresponds to a Type 
lib SN, although of a peculiar kind due to the weakness of the He I 
lines. 

• For the heavily mixed models Bmf4p41x4 and Bmf5p09x4, 
the spectral appearance is markedly different from all other mod- 
els discussed above. In both models, and despite the differences in 
ejecta mass, the spectral appearance is the same and very typical of 
SNe lb in the sense that He I lines are now unambiguously seen at 
all times. Because of the mixing procedure, which affects the en- 
tire ejecta, model Bmf4p41x4 does not show H I lines initially (the 
surface hydrogen mass fraction is too low). If the surface hydrogen 
mass fraction had been kept as high as in model Bmf4p41xl, model 
Bmf4p41x4 would be a Type lib SN, but given the absence of H I 
lines, model Bmf4p41x4 (just like model Bmf5p09x4) corresponds 
to a SN lb. 

So, while all these He-rich ejecta sh ow He I lines at early post- 
explosion times, as discussed in detail bv lDessart et al.l J201 ll) . only 
the moderately and heavily mixed models show optical He I lines 
after about 10 d, throughout the peak of the light curve and beyond. 
This also applies to the He I lines at 10830 A and 20581 A. 

These models suggests that there is a continuum of He I line 
strengths as we increase the mixing efficiency (this would have 
been more obvious if we had had intermediate models between 
Bmf4p41x3 and Bmf4p41x4, for example). As mixing increases, 
non-thermal excitation and ionization become more efficient at and 
above the photosphere (Section|4j. The enhanced ionization popu- 
lates excited levels of He I producing He I lines that are otherwise 
vanishingly weak. In these SN lb models the He I lines dominate 
over weaker overlapping lines in the optical. In contrast, the He I 
lines are weak in our SNe Ic models, and one sees instead a com- 
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Figure 8. Evolution with respect to days since the light curve peak of the radius (top-left), the temperature (top-right), the velocity (bottom-left) and the 
overlying mass at the electron-scattering photosphere for both model sequences Bmf4p41xn and Bmf5p09xn. Notice the non-monotonic photospheric velocity 
in the weakly unmixed Bmf5p09xn models, a feature that disappears as mixing is enhanced. 



plex mixture of weak lines associated with Si II, Fe II, Na I, C II, and 
O I. In our simulations, these weaker lines tend to be always present 
and thus, what yields here the lb or Ic classification depends on how 
strong the He I lines are, which in our set of helium-rich ejecta mod- 
els is more related to the abundance and spatial distribution of 56 Ni 
rather than the He abundance (line identifications for observed SNe 
lb and Ic are given in the next section). 

To illustrate this feature, we show in Fig.QT]the variation at the 
photosphere of the mass fraction of helium and of unstable isotopes 
associated with 56 Ni decay in all models (the time origin is the time 
of light curve peak). Varying the mixing efficiency leads at the light 
curve peak to variations of a few tens of percent for the helium mass 
fraction, but it causes orders of magnitude variations in the 56 Ni 
mass fraction. We find that He I lines only appear strongly if the 
oris inal 56 Ni mass fraction in the corresponding layers is greater 
that ^0.01. We expect some sensitivity to the helium mass fraction 
too but not as dramatic. In the heavily mixed models, a helium mass 
fraction of 0.2-0.3 deep in the ejecta is still enough to produce He I 
lines. 

This 56 Ni mass-fraction threshold makes the production of 
SNe lb a challenge - all ejecta without adequate 56 Ni mixing will 
produce SNe Ic. This also suggests a maximum mass for a SN lb 
ejecta of no more than a few Mq. Indeed, the typical SN Ib/Ic 
56 Ni mass is 0.1 Mq. For a fully mixed 10 Mq ejecta, this would 
produce a 0.01 56 Ni mass fraction. Since even in our fully mixed 



models the 56 Ni mass fraction drops by a factor 5-10 between the 
inner and outer regions, the SN lb ejecta mass cannot be more 
than ~3 Mq for the successful production of He I lines. Because 
56 Ni needs to be quite abundant in those helium-rich layers, SNe 
lb have to be fundamentally associated with lower mass ejecta. As 
we discussed recently, they cannot be associated with high mass 
ejecta because of the short rise-time to light curve peak, the narrow 
light curve peak, and the fast light curve declin e at nebular epochs 
jEnsman & Wooslevll 19881 : iDessart et alfcoill) . 

In the massive-star community, high mass progenitors leading 
to W R stars are often considered as potentia l candidates for SNe 
Ib/c tairschi et ai]|2004 iGeorgy et al]|2009l) , but the correspond- 
ing stellar-evolution simulations yield final masses that are typi- 
cally too large to accommodate the observed SNe Ib/c light curves. 
They are also more generally characterized by an envelope with a 
large oxygen-rich buffer between the core and any remaining outer 
helium-rich shell, potentially compromising the efficiency of mix- 
ing, and hence of non-thermal ionization/excitation. An alternative 
is fast rotation, which can allow lower-mass massive stars to die as 
a moderate-mass WR star with some outer helium-rich shell and 
thus one cannot completely exclude this route. To resolve this is- 
sue requires that all stellar-evolutionary models for such objects be 
evolved until iron-core collapse to permit the modeling of the SN 
light curves and spectra and a direct confrontation to SN observa- 
tions. 
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Figure 9. Left column: Montage of synthetic spectra (black) for the models Bmf4p41xl (left; model Bmf4p41x2 has a similar evolution and is thus not shown), 
Bmf4p41x3 (middle), and Bmf4p41x4 (right) from 2 until 100 d after explosion. For illustration, we overlay the synthetic spectrum that results when excluding 
He I bound-bound transitions in the formal solution of the radiative-transfer equation (red). The left panel correspond to a SN lie evolution, owing to the strong 
Ho line initially, and the transition to a spectrum without lines of H I, He I, nor Si II in the optical after about 15 d. The middle panel corresponds to a SN 
lib evolution, with presence of Ha initially and He I lines, albeit anomalously narrow, at later times. The right panel corresponds to a rather standard SN lb 
evolution, with broad He I lines at all times. 
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Figure 10. Same as Fig. [9] but now for the Bmf5p09xl^l model set. Owing to the larger ejecta mass, only the model with the largest mixing (right) produces 
a genuine SN lb, while the other two series correspond to a SN Ic evolution. 



In helium-rich ejecta, the absence of optical He I lines may 
simply come from a reduction in 5 Ni mass fraction in the helium- 
rich layers, either because of moderate or weak mixing or because 
of an asymmetric o6 Ni distribution that fails to thoroughly mix with 
the helium-rich ejecta layers. At the very least, helium deficiency 
is not a prerequisite for the production of a SN Ic, although it is a 



sufficient condition. Additional explanations may be that some SN 
Ic ejecta are more massive, so that, even if fully mixed, the larger 
mass will cause greater dilution and weaker non-thermal excita- 
tion. Non-thermal processes are not selective - atoms/ions share 
the available non-thermal ionization and excitation energy largely 
in proportion to their relative abundances. 
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To close the discussion on our synthetic spectra, we present 
evidence that mixing must occur at some level in all SNe Ib/c. Rel- 
ative to unmixed models, mixing causes the colors at peak to be 
redder and more c ompatible with the observations (Fig. [12] e.g., 
iModiaz et al .120091) . This issue requires detailed spectral modelling 
since an uncertain reddening can compromise the inferred colors of 
observed SNe Ib/c and synthetic spectral colors are sensitive to the 
completeness of the model atoms. In any case, mixing is likely a 
generic property of Type I CCSN ejecta, and it is probably varia- 
tions in the resulting mass fractions for both helium and 56 Ni that 
control the classification as lb or Ic. 



8 COMPARISON TO OBSERVATIONS 

The goal of this paper is not to do a comprehensive comparison of 
our model results to observations. Rather, we are primarily inter- 
ested in understanding the key physics and ejecta properties that 
influence the light curves and spectral evolution of Type lb and Ic 
SNe. Although we use only two initial ejecta models, it is interest- 
ing to check how well our models fare against some well observed 
SNe Ib/c. 

In Fig. [13] (left panel), we show a comparison of the heavily 
mixed model Bmf5plx4 at 10.3 d after explosion wi th the observa- 
tions of the Type lb SN2008D jModiaz et ai1l2009l) on the 2nd of 
February 2008, which corresponds to 23.88 d after the X-ray detec- 
tion (we choose a time when the model matches both the SN2008D 
spectral color and the line widths approximately). Besides a good 
fit to the He I lines we reproduce well the effects of line blanketing 
in the blue part of the optical, as well as the features due to inter- 
mediate mass elements (IMEs) in the red part of the spectrum. De- 
spite ubiquitous line overlap, He I lines are strong enough to be un- 
ambiguously identified. For example, He 1 5875 A remains stronger 
than the overlapping lines due to Fe II, Si II, and to a weaker extent 
C II. In the near-infrared, He 1 10830 A overlaps with two weak C I 
and O I lines, and a strong Mg II line, but the broader component is 
clearly due to helium. This is in agreement with the chemical strat- 
ification of the ejecta, with a larger abundance of IMEs at smaller 
velocities in such a helium-rich progenitor star. 

In the right panel of Fig. Q~3] we show a similar compari- 
son, but now to th e Type Ic SN 2007gr at ~25 d after discovery 
dValenti et al.l2008l) . using the moderately mixed model Bmf5plx3 
at 15.1 d after explosion. Unlike the spectrum of SN2008D dis- 
cussed above, the He I lines are very weak and dominated by over- 
lapping components. For example, He 1 5875 A is now weak so that 
this spectral region is composed of a number of overlapping lines 
due to Na I, Fe II, Si II, and O I. In the near-IR, the He 1 10830 A line 
is now weaker/narrower and suffers greater contamination from C I 
and Mg II lines. While SN 2007gr is classified a Type Ic, a model 
with weak helium lines (but truly helium-rich), provides a satisfac- 
tory match to observations. This shows that the evidence behind the 
helium deficiency of SNe Ic is weak, if not inexistent. 

A property of SNe lb that is interesting to examine is the 
width of the He 1 5875 A line. For example, for the standard SN 
lb 2008D, the blueshift with respect to line center of the max- 
imum absorption in this line corresponds to a velocity of about 
10000 km s _1 at the light curve peak. In Fig. [14] we show the cor- 
responding velocity versus time since light curve peak for the mod- 
erately and heavily mixed models (we overlay the trajectory of the 
photospheric velocity for each model as a solid line). In the heav- 
ily mixed models, we obtain velocities at maximum absorption in 
He 1 5875 A of about 7000-8000 km s" 1 , which is close but some- 



what short of the observed standard. Interestingly, the moderately 
mixed model Bmf4p41x3, which has the same mass and kinetic 
energy as model Bmf4p41x4, has a maximum-absorption velocity 
in He 1 5875 A about 50% less than that of Bmf4p41x4. This might 
mean that this model is too massive for its energy, or its energy is 
too low for its mass, a combination of both, or that the mixing is 
still too weak. In any case, although not surprising in light of the 
present discussion, this highlights the difficulty of determining the 
kinetic energy of SN lb ejecta, even with the use of spectroscopic 
information. This inference requires the treatment of non-thermal 
effects and detailed radiative-transfer modelling (including time- 
depe ndence since it does als o modify the ejecta ionization struc- 
ture; iDessart & Hillierll2008l) . Such a physically-consistent treat- 
ment is generally not performed, which suggests the inferred char- 
acteristics of SNe Ib/c published in the literature today are not con- 
verged. 



9 DISCUSSION AND CONCLUSIONS 

We have presented results from non-LTE time-dependent radiative- 
transfer simulations based on 1.2 B piston-driven expl osions pro- 



duced from two of the binary-star evolution models of Yoon et al 



20 IC ). This work departs from the previous study of IDessart et al 



20111) by the explicit treatment of non-thermal processes which are 



associated with the high-energy electrons arising from Compton- 
scattering of 7-rays emitted by decaying 56 Ni and 56 Co nuclei. Our 
approach solves for the relative contribution of these non-thermal 
electrons to heating, ionizing, and exciting the ejecta gas, and al- 
lows us to compute accurate spectra and light curves through the 
entire photospheric phase. 

This work primarily focused on exploring the effects of non- 
thermal processes on the photospheric phase of SNe Ilb/Ib/Ic. 
Rather than using a large grid of models with very diverse prop- 
erties, we choose for clarity to focus on only two ejecta models, 
but arising from binary star evolution. For each of these models 
we create four additional models that are characterized exclusively 
by different levels of mixing. This mixing is performed through a 
"boxcar" algorithm which affects all species indiscriminately. We 
arrive at the following conclusions: 

(i) The most important consequences of mixing are associated 
with the change in 56 Ni distribution. We identify both thermal and 
non-thermal effects. 

(ii) Under ionized conditions, and thus generally at optical- 
depths in excess of unity, the non-thermal electrons dump most 
of their energy in the form of heat, thus primarily modulating the 
ejecta temperature, its ionization state, its optical depth, and thus 
the light curve evolution. With enhanced mixing, the light curve 
starts its post-breakout re-brightening phase earlier, follows a more 
gradual rise to peak, and peaks at an earlier time. The post-breakout 
plateau is absent in the most mixed models we study. In all our 
models, 7-rays are trapped for up to ~30d after peak and thus 
ejecta differing only in mixing properties have similar light curve 
slopes in the first month after peak. However, with enhanced mix- 
ing, 7-ray escape starts earlier and causes the nebular luminosity 
to decrease faster. This sensitivity to mixing conditions are exacer- 
bated by the small ejecta masses that characterize all our binary-star 
progenitor models. 

(iii) While non-thermal ionization and excitation are weak in 
optically-thick ionized layers, non-thermal electrons in partially- 
ionized regions dump a significant fraction of their energy through 
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Figure 11. Evolution of the helium (black), oxygen (blue), and 56 Ni (including its decay products; red) mass fractions at the photosphere for the simulations 
associated with the Bmf4p41xl^t (top) and Bmf5p09xl^l (bottom) model sets and shown with respect to the time of peak bolometric luminosity. All 
"standard" SNe lb show strong He I lines at maximum brightness. Only models Bmf'4p41x4 and Bmf5p09x4 have this property, which seems to require a 56 Ni 
mass fraction on the order of 0.01 in the helium-rich regions. 
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Figure 12. Spectral comparison from UV to near-IR and at light-curve peak for models Bmf4p41xl— 4. A distance of lOMpc is adopted. With increasing 
mixing, colors are redder (the photosphere is cooler and more strongly affected by line blanketing in the UV), line profiles are broader (the strongest lines are 
optically thick to larger distances/velocities), and the photospheric velocity is larger (the continuum also forms further out). 



non-thermal ionization and excitation of the plasma. In heav- 
ily mixed models, the significant decay energy deposited above 
the photosphere affects primarily the He I ionization and excita- 
tion since helium is nearly recombined there and is the dominant 
species. In contrast, weakly or moderately mixed models deposit 
most of the energy at depth, below the photosphere, and fail to non- 
thermally excite (or ionize) the outer ejecta region where lines and 
continuum form. However, in these models, the exclusive deposi- 
tion of heat at large optical depths creates highly ionized condi- 
tions. 

(iv) Consequently, during the photospheric-phase of SNe Ib/c, 
no n-thermal eff ects are favored in highly mixed ejecta. As proposed 
by ILucvI dl99ll) . non-thermal excitation and ionization are key for 
the production of He I lines. In our present set of helium-rich mod- 
els, the primary feature is the presence of optical and near-IR He I 
lines characteristic of the Type lb classification. Although charac- 
terized by the same cumulative yields, moderately mixed models 
fail to non-thermally excite helium atoms and lead to a Type Ic 
classification. 

(v) In spite of our small model set, we identify a diverse range 
of spectral evolution and classification, and obtain both Type lb and 
Ic SNe. Furthermore, some models show hydrogen at e arly times 
(produced by non-LTE effects alone; iDessart et al.l201 ll) . followed 
by the presence or absence of He I lines. These would correspond 
to classifications as Type lib and lie, resp ectively. The latter has 
never been identified with confidence (see iBranch et al.ll2006l for 
a possible identification of Ha in the Type Ic SN 19941) but this 



may just be a bias of observational campaigns which rarely capture 
Type Ib/c at a few days past explosions. In view of the present sim- 
ulations, a type lie classification is easily understood as stemming 
from a progenitor with both hydrogen and helium, the former be- 
ing under-abundant but easily seen, and the latter arbitrarily abun- 
dant but quickly inv isible in the absence of non-thermal excitation 
jDessart et alj|201 lb . A transition from type Ic to type lb is more 
difficult to achieve as it would require a large amount of helium 
together with moderate mixing. Non-thermal excitation of helium 
may then become effective only when the photosphere has receded 
to deeper layers where 56 Ni is mor e abundant. This uniq ue transi- 



tion has been seen for SN2005bf bv lFolatelli et~aU(l2006l), althoug h 



— ^ 

mixing is unlikely the cause in this special case l IMaeda et al.l2007l) . 

(vi) Apart from the cases where He I lines are really strong, line 
overlap considerably complicates the identification of spectral fea- 
tures. We identify a slew of lines from He I, C I, O I, Na I, Mg II, and 
Fe II in the optical and the near-IR. One can, for example, easily 
mistake He 1 5875 A with Na I D, Si II 6350 A with Ha etc. Multi- 
epoch spectroscopic observations from a few days after explosion 
until the nebular phase are critical to build a comprehensive and 
accurate picture of the SN-ejecta composition. 

(vii) From our explorations, we are led to associate SN lb ejecta 
with the explosion of low-mass helium cores in which efficient 
small-scale mixing takes place. The scale of relevance for non- 
thermal excitation/ionization is the 7-ray mean-free-path, which 
varies with both depth and time. We favor progenitors that have 
only a small oxygen-rich mass buffer between the 56 Ni produc- 
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Figure 13. Left: Comparison of model Bmf'5plx4 at 10.3 d after explosion (we show the total synthetic flux in red and the synthetic continuum flux in blue), 
redde ned by E(B — V ) = 0.7, with the observations of the Type lb SN2008D on the 2nd of February 2008, 23.88 d after the X-ray detection l lModiaz et all 
|2009|) . Right: Similar as the left figure, but now showing a comparison of model Bmf5plx3 at 15 .1 d after explosion , reddened by E(B — V) = 0.15, with 
the observations of the Type Ic SN2007gr on the 10th of September 2007, ~25 d after discovery dValenti et all2008h . Here, we illustrate the correspondence 
between these observations and our models when both are characterized by the same color and line widths. In each panel, the lower part shows the relative 
contribution to the total spectrum of bound-bound transitions associated with individual species. Notice in particular the distinct He I contribution in the 
strongly mixed (left) and the moderately mixed model (right). 



tion site and the outer helium-rich layers, facilitating the mixing, 
and preventing too much dilution of 5f> Ni into a large mass. In our 
set of models, we find that a minimum 5 Ni mass fraction of 0.01 
is necessary at large velocities to trigger significant non-thermal 
effects early on. The small-scale mixing necessary in this con- 
text may originate fro m the fragmented structure of the SN shock 
dNordhaus et alj2010h . 

(viii) It requires special circumstances for our low-mass quasi- 
hydrogen-less cores to become a lb — failure to meet these require- 
ments will typically lead to a type Ic. This may stem from any com- 
bination of properties that prevent helium atoms from being non- 
thermally excited. This difficulty to make SNe lb might be reflected 
by the 2:1 occurrence rates of SNe Ic:Ib dSmarttll2009l) . Inefficient 
non-thermal excitati on of He I may resu lt from 1) helium deficiency 
or under-abundance dYoon et alJCOlol argues for the possibility of 
even lower mass helium-deficient progenitors than those presented 
here); 2) a higher-mass helium core with a sizeable oxygen-rich 
mass buffer, which causes a stronger dilution of 56 Ni — this could 
naturally occur for single WR stars which are born from high- mass 
stars with higher-mass helium cores dGeorgy et al.l2009-l) ; 3) explo- 
sions that are not mixed on small scales so that most of the 56 Ni 



is contained within a limited solid-angle and ejecta volume. In this 
context, a "jet" of 56 Ni-rich material is not expected to be a suitable 
morphology for non-thermal excitation of ejecta material because 
the 7-rays would not be able to penetrate the other latitudes that 
would contain the bulk of the helium-rich material, i.e., only a tiny 
fraction of He I atoms would be excited. 

(ix) Because the combination of thermal and non-thermal pro- 
cesses considerably affect ion/atom excitation and thus the asso- 
ciated line optical depths, we find that models of the same mass 
and kinetic energy but with enhanced mixing display broader line- 
profile widths at a given post-explosion time. This sensitivity com- 
plicates the reliable inference of SN-Ib/c ejecta masses and ex- 
plosion energies, and suggests the importance of accounting for 
both time-dependent and non-thermal effects when modelling their 
spectra and light curves. 

(x) Adding to co nclusions of earl i er studies 

dEnsman & Wooslevl 1 19881 ; IShigevama etal] Il990l : | Dessart et al.l 
1201 lh . our spectral simulations including non-thermal effects 
strengthen the notion that the bulk of SNe Ib/c arise from low-mass 
ejecta likely resulting from binary-star evolution. The broad line 
Ic's with and without GRBs, associated with larger mass ejecta 
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Figure 14. Time evolution of the velocity at maximum absorption in the 
strongest optical He I line at 5875A. We only select here the models that 
show weak (Bmf4p41x3) and strong (Bmf5p09x4 and Bmf5p09x4) He I 
lines. We also overlay in each case the evolution of the velocity at the 
electron-scattering photosphere (solid line). 



and higher kinetic energy dWooslev et al.l IT999b . may stem from 
different progenitors and different explosion mechanisms. Our 
simulations suggest that these may in fact contain some helium, 
which is invisible because it is not excited. This result does not 
conflict with the predi ctions of massive s tar evolution for carbon 
and oxygen WR stars dGeorgv et al ] |2009h . Furthermore, it would 
not co nflict with the prerequisites of the colla psar model ( Woosleyl 
19931) or those of th e proto-magnetar model dWheeler et al.l l2000; 
Metzger et ap201ll) of long-duration GRBs. 

(xi) We identify a negative feedback on non-thermal effects. As 
the 56 Ni abundance is increased, the heating caused by energy de- 
position increases and forces non-thermal electrons to channel their 
energy in the form of heat, thereby quenching the non-thermal ion- 
ization and excitation channels. In contrast, non-thermal effects 
thrive under neutral conditions, i.e., plasmas that are cool, and 
thus not too 56 Ni-rich. We thus anticipate non-thermal ionization 
and excitation to be stronger for moderate 56 Ni abundances. They 
should be more efficient in the cool photospheres of SNe Ib/c than 
in the hot ionized photospheres of SNe la, highly 56 Ni-rich, where 
non-thermal electrons will tend to dissipate their energy as heat. 

(xii) Our radiative-transfer modeling is ID so the emergent spe- 
cific intensity or flux is independent of viewing angle. As we inde- 
pendently argue in Sect. [2] multi-dimensional effects are unlikely 
to turn a SN lb into a SN Ic, or a lib into a lb, but instead should in- 
troduce slight variations in line width and/or strength with viewing 
angle, as illustrated in o ur ID models dif fering in mixing magni- 
tude (Fig 1 14b. Recently, iRest et al] 1201 lh have used light echoes 
to reveal the viewing-angle dependent intensity from the type lib 
SN associated with Cas A, typically characterized by variations in 
line widths along different directions. However, each recorded light 
echo sup ports the lib classification for the SN. Using a similar ap- 
proach to Dess art & Hillietl ( 1201 lbl) . we will investigate the magni- 
tude of the asymmetries required to explain such observations. 

The next step in our modeling effort is to compute light curves 
and spectra for a large grid of progenitor stars, both single and from 
binary systems, exploded with a range of explosion energy, and 
characterized by different levels of mixing. Such a diversity of ini- 
tial conditions should reflect the diversity of SNe Ilb/Ib/Ic seen in 



Nature and allow a better assessment of the progenitor and explo- 
sion properties of these CCSNe. 
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